A COMPARISON OF METAL ENRICHMENT HISTORIES IN RICH 
CLUSTERS AND INDIVIDUAL LUMINOUS ELLIPTICAL GALAXIES 



Fabrizio Brighenti 2 ' 3 and William G. Mathews 2 
2 University of California Observatories/Lick Observatory, Board of Studies in Astronomy and 
Astrophysics, University of California, Santa Cruz, CA 95064 
mathews@lick.ucsc.edu 
3 Dipartimento di Astronomia, Universita di Bologna, via Zamboni 33, Bologna 40126, Italy 

brighenti@astbo3.bo.astro.it 

ABSTRACT 

The large spatial extent of hot, X-ray emitting gaseous halos around massive elliptical 
galaxies indicates that most of this gas has not been generated by stellar mass loss. Instead, 
much of this gas results from an intergalactic gaseous inflow toward the overdensity from which 
giant ellipticals and their associated galaxy groups formed. Since these hot gaseous halos are 
old, they contain important information about the star formation history of elliptical galaxies. 
In this paper we show that the enrichment history of this hot gas is closely linked to its gas 
dynamical history; supernovae provide both energy and metal enrichment. We find that Type II 
supernovae based on a Salpeter IMF, plus a small number of additional Type la supernovae, can 
explain the the density temperature and abundance profiles currently observed in gaseous halos 
around massive ellipticals. Within the central, optically bright region of luminous ellipticals, 
approximately half of the interstellar iron is produced by Type la supernovae and half by mass 
lost from evolving stars which were originally enriched in iron by Type II supernovae. However, 
iron and silicon abundances in the intracluster gas within rich clusters suggest enrichment by 
a larger number of supernovae per unit optical light than we require for massive ellipticals. 
The additional supernovae implied by cluster data cannot be reconciled with our models 
for individual massive ellipticals. Evidently, rich clusters cannot be constructed by simply 
combining ellipticals and their associated groups since the enrichment history of clusters and 
massive ellipticals appears to be fundamentally different. Neither currently discussed resolution 
of this discrepancy - increased number of Type II supernovae (flat IMF) or strong Type la 
enrichment in clusters - is attractive. Although the global hot gas iron abundance is similar 
in all large galaxy clusters, silicon is enhanced in hotter, richer clusters. This Si/Fe variation 
implies that E and SO galaxies are not the only sources of cluster gas enrichment; perhaps spirals 
or low mass galaxies are also important. 



Subject headings: galaxies: elliptical and lenticular - galaxies: formation - galaxies: evolution 
galaxies: cooling flows - x-rays: galaxies 



INTRODUCTION 



Luminous elliptical galaxies are the astronomical analog of tree rings or geological core samples - within 
the hot gas in these massive galaxies resides important information about the earliest stages of galaxy 
formation and the enrichment of the cosmos with heavy elements from the first generations of stars. The 
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stellar density profile and low rotation of bright ellipticals indicate that they were assembled from smaller 
galactic objects containing stars. The most likely sites for the formation of giant elliptical galaxies are 
small groups of galaxies where mutual tidal interactions are strong and the likelihood of galactic merging 
is enhanced. In an idealized version of this formation hypothesis, the merging process within the group 
comes to completion after a few dynamical times when the massive central elliptical has consumed the dark 
and luminous matter of many original group galaxies of moderate mass. When this happens, the merging 
process shuts down and the elliptical evolves passively as a group-dominant elliptical in the field or enters 
a large cluster. While there is clear evidence that some massive ellipticals contain younger stars and must 
have suffered significant mergers since the epoch of most intense galaxy formation, there is also evidence 
that other more venerable ellipticals have survived until the present time, relatively undisturbed by events 
that have occurred since their formation at early times. The hot interstellar gas in these old galaxies is 
most interesting since it still contains information about the earliest stages of star formation, including the 
heating and enrichment of the gas by supernovae. 

In this paper we investigate the production of iron and silicon by Type la and Type II supernovae 
(SNIa and SNII respectively) and follow the dynamical redistribution of these metals in the hot interstellar 
and intragroup gas until the present time. Gas enrichment and dynamical evolution are linked since the 
collective energy released by Type II supernovae can expel metal-enriched gas from elliptical galaxies in a 
wind. The energetics of Type II supernovae during early star formation can be constrained by comparing 
heavy metal abundances expected from evolutionary gas dynamical calculations with abundances provided 
by X-ray observations of individual massive ellipticals. In a similar fashion, metal abundances observed in 
rich clusters of galaxies provide global enrichment constraints on early type galaxies which are generally 
regarded as the principal source of metals in the intracluster gas. Rich clusters are thought to have retained 
all of the products of supernova enrichment either within galactic stars or in the hot intracluster medium. 
In this sense rich clusters are closed boxes. From an observational standpoint individual massive elliptical 
galaxies are certainly not closed boxes since most of the metals produced by supernovae have been carried 
out by supernova-driven winds. Much of these metals now resides in low density gas at large distances 
from the optical galaxy where detection is difficult or where the metal enriched gas has been tidally or 
dynamically removed. From the standpoint of our gas dynamical models, however, individual ellipticals are 
closed boxes since we can accurately follow the products of supernova enrichment into distant regions of 
very low density. 

Thermal X-ray emission from hot circumgalactic gas surrounding many bright ellipticals often extends 
far beyond the optical image of the galaxy (Mathews & Brighenti 1998a). The large mass and extended 
spatial distribution of this gas cannot be understood only from normal evolutionary stellar mass loss. Our 
dynamical models indicate that the most distant gas in these galactic halos has accumulated by secondary 
infall into the perturbation that initially led to the formation of the group and its dominant elliptical 
(Brighenti & Mathews 1998a, 1998b; Mathews & Brighenti 1998b). 

The objective of our recent gas dynamical calculations has been to follow the evolution of the hot 
gas within and around massive elliptical galaxies since the time of galaxy formation. A successful model 
can reproduce the radial density and temperature distributions as determined from X-ray observations. In 
our recent models we begin with an overdensity perturbation in a simple flat cosmology. As dark matter 
accumulates in the core of the growing perturbation, we assume that a stationary NFW dark matter halo is 
formed (Navarro, Frenk and White 1996). Exterior to this stationary halo the dark matter inflow is identical 
to the cosmic similarity flow of Bcrtschingcr (1985). The baryonic gas evolves in this time-dependent 
potential: flowing inward, shocking at some radius and radiatively cooling near the center. At some early 
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time, t^, when enough baryons have accumulated in the perturbation, we form the de Vaucouleurs profile 
of the large elliptical observed today. By this means we circumvent the complex merging processes that 
occurred at very early times. Nevertheless, our calculation globally conserves dark and baryonic mass and 
treats the gas dynamics in full detail from t. t to the present time. In these simple models we release the 
Type II supernova energy at the same time £» that the stellar system is formed but the Type la energy is 
released from £* to the present according to some assumed variable rate. 

However, there is a degeneracy in our successful solutions. We find that the hot gas density and 
temperature profiles typical of massive ellipticals, n(r) and T(r), can be reproduced with different 
combinations of the fundamental model parameters: the time of star formation £#, the total SNII energy 
released En, and the universal baryon mass fraction f2{, « 0.05 ± 0.01 as determined by big bang 
nucleosynthesis (Walker et al. 1991). For example, solutions with 0,^ decreased by <~ 0.01 are similar to 
those with En increased by ~ 2 or with i„ decreased by ~ 1 Gyr. 

Much of the degeneracy involving these uncertain parameters can in principle be removed by considering 
the abundance and of iron and other elements produced in supernovae. The total mass of heavy elements in 
stars and gas today is directly related to the total number of supernovae of each type that have occurred. 
Unfortunately, many of the essential supernova-related parameters are poorly known: the average amount 
of iron or silicon created in each supernova event, the total mass of these elements locked into stars today, 
the initial mass function (IMF), etc. 

While we cannot attempt here a definitive theoretical resolution of the enrichment history of early type 
galaxies, we can make progress within the limitations of current theory and observational data. In view 
of the uncertainties involved, however, a highly detailed gas dynamical calculation including the effects of 
metal enrichment seems unwarranted. The approximate results we present here can be regarded as the first 
step in an iterative procedure that will become better defined in the future. 

Nevertheless, the metal enrichment history of elliptical galaxies that we describe here is sufficient 
to rule out several scenarios that have been widely discussed. One of our motivations for beginning this 
project is the apparent dissimilarity of iron abundance and gas fraction between galaxy groups and rich 
clusters which led David (1997) and Renzini (1997) to remark that present-day, rich clusters cannot be 
assembled from present-day groups. For example, most (50% - 80%) of the mass of heavy metals in rich 
clusters resides in the intracluster gas where the iron mass per unit of stellar light Mf s /Lb consistently 
exceeds that in small groups or individual massive ellipticals. Moreover, the lower fraction of gas mass 
in galaxy groups relative to clusters indicates that much of the metal-enriched gas produced in groups or 
individual ellipticals has been expelled into the environment by the release of supernova energy (Arnoud, 
et al. 1992; David, Jones, & Forman 1995). Perhaps, therefore, the gas and iron formerly within large 
ellipticals and galaxy groups has been expelled by supernova-driven winds into the low density environment 
where it escapes X-ray detection. However, this low density, metal-enriched gas can be followed in our 
gas-dynamical models. In successful models the metal rich gas remaining within the galaxy must agree with 
the gas abundance and abundance gradients observed in these bright galaxies. 

Our results reported here indicate, however, that the mass of metals that has flowed out of luminous 
ellipticals into their low density environment is insufficient to account for the larger mass of metals observed 
in rich clusters. We are therefore faced with a fundamental inconsistency: parameters that allow the 
intracluster gas to be enriched by galactic winds from massive ellipticals are inconsistent with detailed 
X-ray observations of the metal content within individual ellipticals. 

The relative importance of SNII and SNIa in creating the metals observed in cluster gas has been 
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widely debated. The total iron abundance in rich clusters is too large to be produced by SNII with a 
standard Salpeter IMF. Many authors have suggested that SNIa produce enough iron to account for the 
additional iron required (e.g. Ishimaru & Arimoto 1997; Renzini 1997; Wyse 1997). This solution is 
appealing since approximately | of the iron in our Galaxy is produced by Type la supernovae (Renzini 
1997). However, we show here that the total silicon produced by both types of supernova, when considered 
along with additional gas dynamical constraints, indicates that SNIa produce at most only a few percent of 
the total iron in massive ellipticals. By contrast, if SNII are the primary source of metal enrichment in rich 
clusters then the IMF must be flatter than Salpeter; this solution to the iron problem has been suggested 
by David (1997) and Gibson, Loewenstein & Mushotzky (1997). However, the iron abundance and its 
radial gradient observed in massive ellipticals cannot be fit with our gas-dynamical models if the number of 
SNII exceeds the Salpeter value. It would appear therefore that massive ellipticals may not be the primary 
source of enrichment of the intracluster gas. 

We begin our discussion with a review and reanalysis of the observed iron and silicon abundances in 
rich clusters and discuss the implications of these abundances for clusters and individual early type galaxies. 

2. SUPERNOVA ENRICHMENT IN RICH CLUSTERS 

The mean iron abundance observed in the hot gas in rich clusters is about 0.40 in solar units with 
a cosmic scatter (0.3 - 0.65) that exceeds the errors of measurement (Mushotzky 1998). These values 
are somewhat greater than the iron abundances quoted in the X-ray literature, most of which are based 
on the "photospheric" solar iron abundance (by mass) z^ P p\ Q = 2.66 x 10~ 3 since we assume here the 
lower, "meteoritic" , solar iron abundance ZFe,o = 1-83 x 10~ 3 . This distinction between "meteoritic" and 
"photospheric" solar iron abundance was made in the review by Anders & Grevesse (1989). However, recent 
atmospheric models of the sun indicate that the iron abundance in the solar photosphere is similar to the 
so-called "meteoritic" value of Anders & Grevesse (Holweger et al. 1991; Biemont et al. 1991; McWilliam 
1997). The implications of the lower zf s ,q for X-ray abundances has been discussed by Ishimaru & Arimoto 
(1997) and Gibson, Loewenstein, & Mushotzky (1997). 

Arnaud et al. (1992) have shown that the total gas mass in the hot intracluster medium (ICM) of rich 
clusters A4 g and the mass of iron in that gas are proportional to the total optical luminosity Cb of all early 
type galaxies in the cluster. Arnaud et al. find that M g is uncorrelated with the total Cb,sp from spiral 
galaxies in these clusters. Based on these findings, it is generally assumed that early type galaxies dominate 
the cluster gas enrichment. 

The relative contributions of Type II and Type la supernovae (SNII and SNIa) to the ICM metal 
enrichment has been widely discussed. The total mass of a particular heavy element in the ICM can be 
determined from the mass of that element produced by supernovae and the total number of supernova 
events in the past. The number of SNII events can be inferred from the total spheroidal stellar mass in 
the cluster, together with an assumed initial mass function (IMF) to determine the fraction of stars that 
become Type II supernovae. For SNIa the currently observed rate for ellipticals is poorly known and 
almost nothing is known about the SNIa rate at earlier times. The SNIa rate is thought to depend on mass 
exchange between binary stars and theory provides no unambiguous way to estimate orbital parameters 
and binary mass ratios. The yields of heavy metals in SNIa are thought to be known reasonably well 
from theoretical models while those from SNII are uncertain because of complexities in the pre-explosion 
evolution, uncertain reaction rates, and the unknown amount of processed material that falls back on the 
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stellar remnant (Woosley & Weaver 1995; Gibson, Loewenstein & Mushotzky 1997). 

In view of these difficulties, the relative influence of SNII and SNIa in enriching elliptical galaxies and 
rich clusters must be determined largely from observations of the metal abundance in stars and gas. The 
enrichment history of individual ellipticals and associated groups of galaxies is difficult to access since the 
total energy released in all SNII is sufficient to expel a significant fraction of metal-rich gas from these 
systems. However, rich clusters of galaxies are massive enough to have retained all of the material processed 
in supernovae over cosmic time. The "closed box" nature of rich clusters is supported by the relative 
constancy of the ratio of baryonic to total mass for cluster masses > 5 x 10 13 M© (David 1997). The 
data assembled by Renzini (1997) indicates that metal enrichment by supernovae in rich clusters remains 
confined within the cluster potential; M g / Cb and the total iron mass per unit optical light MFe,g/ £-b are 
nearly constant for cluster masses > 5 x 10 13 M©. Clusters and galaxy groups with masses less than this 
have lower gas to stellar mass ratios, suggesting gaseous outflow. 

In determining the past history of supernova enrichment, the silicon abundance is more useful than the 
iron abundance. The silicon abundance has been observed in the ICM of many rich clusters and is based 
on a well-understood K-line emission feature so that translation into abundances should be reliable. In 
addition, silicon yields from SNII and possibly SNIa may be better determined from theoretical supernova 
models than those of iron. While the iron yield from SNII yFe,u{m) increases monotonically with the 
pre-supernova stellar mass m, the silicon yield ysi.u{m) has a pronounced maximum at m <~ 22 M© 
provided the energy released by the supernova is not too large (A and B models of Woosley & Weaver 1995). 
Because of this maximum, the IMF-averaged yield (ysiji) is very insensitive to the choice of parameters 
for power law IMFs provided the upper mass limit m u > 22. 

Assuming that silicon is a reliable tracer of past supernova activity in rich clusters, we wish to construct 
an expression for the observed ratio of the total mass of silicon in the hot cluster gas (ICM) to the total 
optical luminosity of spheroidal system stars in the cluster, 

= Msi,g/C B - 

The present mass of silicon in the ICM is determined by the total mass of silicon produced by SNII and 
SNIa less the amount of silicon currently within stars: 

™ _ Mst.II + Msija - M Sl .* n , 
' 9 = ~L~b ~' 

Another important observed quantity is the currently observed silicon to iron ratio in the ICM, 

R=( ZSi \ = ( M Si \ = M Si,II + M S Ua - M Sl ,* ,^ 
~\ZFe) g \Mfc) g MfcJI + M F e,Ia - M F e.*' 

Note that R is the ratio of silicon to iron masses in absolute units and is not normalized with the solar 
ratio. We now seek expressions for Tsi.g and R in terms of the expected number of supernovae and the 
average silicon and iron yields for each supernova event. 

If we assume that all stars more massive than 8 M Q produce SNII, the number of SNII per M© of 
stars formed, rju, can be found for any assumed IMF. For a power law IMF <j)(m)dm — <j) m~( 1+x ' > dm the 
number of SNII is equal to the total number of stars more massive than mg = 8 M©, therefore 

N H x - I mj 1 - m~ x 
mi = M7 = x m]- x - m\- x ' 
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For example, r\ u = 6.81 x 1(T 3 SNII per M© for a Salpeter IMF (slope x = 1.35) having upper and lower 
masses of m u — 100 M© and mi — 0.08 M© respectively. In recent discussions of the metal enrichment of 
early type galaxies and galaxy clusters it has been fashionable to consider only single power law IMFs (e.g. 
Loewenstein & Mushotzky 1996; Ishimaru & Arimoto 1997; Gibson, Loewenstein & Mushotzky 1997) with 
the Salpeter slope regarded as "normal." However, many lines of evidence in our own Galaxy suggest that 
the normal IMF has some curvature, becoming natter at subsolar masses (Leitherer 1998). Scalo (1998) 
suggests the following triple power law approximation for a universal IMF: 



Evaluating the specific SNII frequency for this IMF we find j] n = 7.81 x 10~ 3 which differs by less than 
15 percent from the single-slope Salpeter value. In view of the many larger uncertainties involved in other 
parameters, we shall continue to use simple power law IMFs in this paper to determine the number of SNII, 
mean Type II supernova yields and the stellar mass loss rate. 

The total number of SNII explosions during the entire history of rich cluster stars is A/jj = f]n A4* 
where A4* is the total initial stellar mass in cluster E + SO galaxies. The current mass of all early 
type galaxies in a rich cluster can be estimated from the total B-band luminosity of all cluster galaxies, 
M* = (M/Lb)Cb, where M/Lb ~ 7 is an average mass to light ratio for bright cluster ellipticals. 
Assuming that star formation in ellipticals can be approximated with a single burst at time i#, the stellar 
mass at that time was A4* — A4*/(l — (3) where (3 is the fraction of the initial mass that has been lost 
from the stars between t* and the present time t n . For example, if stars form with a Salpeter IMF at time 
t* = 2 Gyrs and t n = 13 Gyrs is the present time, then (3 = 0.3. 

Since galactic stars in rich clusters are observed to be enriched by SNII ejecta, it is clear that the single 
burst approximation is only an approximation. However, if most of the stars are formed in several bursts 
concentrated near time t#, the long term stellar evolution will be little altered although the single-burst 
value of (3 will be overestimated if some supernova-processed gas is formed into stars in nearly simultaneous 
multi-bursts. If a fraction /* of SNII- produced metals formed into stars at time t*, then a fraction 

= /*(1 — (3) is still locked in stars today. The total mass of SNII-produced silicon still remaining in the 
intracluster gas phase today is therefore 



where the mean silicon yield (ysi.u) is the IMF-averaged mass of silicon (in solar masses) generated per 
SNII event. The total mass of iron in the ICM from SNII is given by the same expression with (ysiji) 
replaced with (yFeji)- By using the single burst assumption, modified by allowing some of the SNII ejecta 
to form stars, we are neglecting higher order details such as second generation SNII formed from stellar 
ejecta. 

The total amount of silicon produced in SNIa is the product of the total number of Type la supernovae 
and the silicon yield, Afi a ysi,ia- Since the past rate of SNIa explosions in ellipticals is poorly known, we 
represent our ignorance with a power law: 



4>(m)dm = fam ( - 1+x ^dm i = 1, 2, 3 



with 



x x = 0.2 ± 0.3 for O.K to < 1 M© 
xi = 1.7 ± 0.5 for 1.0 < to < 10 M Q 
x 3 = 1.3 ± 0.5 for 10 < to < 100 M© 



(i - p)~ 1 (M/L B )CBmi(ys l ji)(i - f,) 



dATia 
dt 



10 10 lOOyrs 



SNu(i) SNIa yr 



,-i 
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whcrc Cb is in units of Lb,q and 

SNu(i) - SNu(t n ){t/t n )- p 

is the SNIa rate in SNu units (supernovae per 10 10 L B per 100 years). The total number of SNIa since t. t is 
therefore 

■A/ja = ^ 1^01^ "/ (t/t n )- p d(t/t n ) = C B n Ia 

where n/ a is the number of SNIa per unit L b ■ The important parameter that governs the SNIa contribution 
to the observed ratios R and Ysi, g is nj a although we shall occasionally use the current SNIa rate SNu(t n ) 
to characterize n/ , assuming particular values of p and t*. If a fraction of the enriched SNIa ejecta was 
incorporated into stars at time i*, then a fraction G* = — /?) remains in stars at the present time. 

Finally, the current stellar iron abundance is given by 

, Mfc* ,„,. 

where ZFe,* is the abundance of iron by mass in stars relative to hydrogen and ZF e .*j\A is the ratio of iron 
mass to total mass including helium. 

Following the notation introduced above, Equations (1) - (3) can be written as 

T S i, g = = (1 - l3)- 1 (M/L B )f lH (y Sl .ii)(l - F t ) + n Ia y SUa (l - G*), (4) 

L-B 

R _ (M Si \ (1 - l3r 1 (M/L B ) mi (y Sl .ii)(l - F*) + n IaVSlJa {l - G*) 

\M F e) q (l-/3)- 1 W^)'7//(yF e ,//)(l-^)+n/ a yFe,/a(l-G*)' U 



and 



ZFe 



1.4[(1 - Py'viliVFeJ^F* + {M/LB^maVF^aG*]. (6) 



If SNIa do not contribute to the enrichment of the stars, then the ratio of stellar silicon and iron must be 
in proportion to SNII yields, z Fe ,* = z Si,*{yFe,n) I '{ysi,ii)\ equation (6) reduces to this simple expression 
when nj a = or G» = 0. With the total number of SNIa regarded as a parameter, we shall solve the three 
equations above for the most uncertain remaining parameters: the specific number of SNII, rjjj, the mean 
iron yield from SNII (yFeji), and the fraction /* of SNII ejecta that has formed into stars. 

Values for the many additional parameters in the last three equations must be determined from 
theoretical supernova calculations or from observation. IMF-averaged theoretical supernova yields for iron 
and silicon have been collected and discussed by Loewenstein & Mushotzky (1996) and Gibson, Loewenstein 
& Mushotzky (1997). Following these authors we take SNIa silicon and iron yields from Model W7 of 
Thielemann, Nomoto & Hashimoto (1993): ysi,i a — 0.158 and yFeja = 0.744, both in M Q . Yields for 
SNII collected by Loewenstein & Mushotzky (1996), mostly based on Woosley & Weaver (1995), have a 
large range for different models: (ysi.n) = 0.08 — 0.32 and (yFeji) = 0.11 — 0.34 while yields for the 
less extreme SNII models discussed by Gibson, Loewenstein & Mushotzky (1997) have a more limited 
range: (ysi,n) = 0.104 - 0.143 and (yFe.ii) = 0.073 - 0.141 both evaluated with a Salpeter IMF. Gibson, 
Loewenstein & Mushotzky (1997) discuss the uncertainties in theoretical SNII yields in detail. As discussed 
previously, because of the maximum in the silicon yield near progenitor mass 22 Mq, the mean silicon yield 
(ysiji) ~ 0.133 Mq appears to be more securely known than (yFeji) which we regard as an unknown to 
be determined by solving the equations above. 
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The stellar mass to light ratio in early type galaxies is a slowly increasing function of galactic 
luminosity, M/Lb oc Lb°' 2 (Faber et al. 1984), and depends somewhat on the cluster luminosity function. 
For a representative value we choose M/Lb = 7 (in solar units), characteristic of typical bright ellipticals. 

Values of Ts i;£ , and R in Equations (4) and (5) must be supplied from X-ray observations. For three 
rich clusters (A 2199, A 496, and AWM 7) observed with ASCA Mushotzky et al. (1996) find an average 
silicon mass to light ratio Tsi.g — 0.0305 ± 0.008 in solar units. More extensive observational data indicates 
that there is a significant cosmic variation in the abundances of Si and Fe among rich clusters (Mushotzky 
& Loewenstein 1997; Mushotzky 1998). If real, this variation of galaxy-averaged values suggests that the 
IMF, star formation efficiency, supernova frequency or other supernova properties may vary dramatically 
among cluster galaxies. Nevertheless, for our purposes here we simply use average abundances from all of 
the clusters observed. ASCA observations indicate a silicon to iron ratio of Si/Fe = 2.2 ± 0.25 (Mushotzky 
1998) in solar units. But this abundance ratio, based on the so-called "photospheric" solar iron abundance, 
becomes Si/Fc = 2.2(1.83/2.66) = 1.51 in units of the "meteoritic" solar ratio (Si/Fc) Q = 0.550 adopted 
here. The corresponding absolute Si/Fe abundance ratio in equation (5) is therefore R = 0.83 

Solutions of Equations (4) - (6) for 77//, {yFe.ii), and /* are shown in Figure 1 as functions of the total 
number of SNIa events A/} and the current Type la supernova rate SNu(i„) evaluated with p = 1, t n = 13 
Gyrs, and i* = 2 Gyrs. Two representative solutions are illustrated: (i) with no enrichment of stars by 
SNIa ejecta g* = and (ii) with a significant SNIa enrichment = 0.5. The plotted solutions are based on 
a mass return parameter appropriate for a single burst Salpcter IMF (3 = 0.3. 

Type la supernova rates observed in elliptical galaxies favor low values of SNu(t„) where ijjj, /* 
and (yFeji) in Figure 1 are almost independent of SNu(t„) and A/} a <C A/jj. In a recent study of 
the observed frequency of SNIa in ellipticals, Cappellaro et al. (1997) find a low current SNIa rate 
SNu(t„) = 0.058(if/50) 2 in SNu units. The current SNIa rate is also restricted by our gas dynamical 
models for the evolution of hot gas in large elliptical galaxies. The computed radial variation of gas density 
and temperature agree with those observed only when SNu(i„) ^ 0.25 (with p = 1), again suggesting low 
values for A// a - However, when silicon and iron abundances are included in these hydrodynamical evolution, 
as discussed below, low values of SNu(i„) and A/} a are essential. In this low-SNu(£„) range of Figure 1, 
solutions for the three parameters r/jj, (yFeji), and /* approach the limit of no SNIa enrichment, A// a — * 0. 
In this limit solutions to Equations (4) - (6) for Tsi, g , R, and z»Fe simplify to 

(VFeJl) = (VSi,Il)/R = 0.16, (7) 

(M/L B )(y Si ji){l-l3)- 1 1 ; 

and 

/. = ^T^? = (9) 

vii\ysi,ii) 

These numerical evaluations are based on Tsi, g — 0.0305, R = 0.83, zp e ,* = 0-§ZFe,Q, M/Lb = 7, 
(ysiji) = 0.133, and (3 = 0.3. Since Ts i ff dominates the numerator of Equation (8), the value of rju is 
essentially unchanged if the mean stellar iron abundance is taken to be solar, z Fe ,* — z Fe,Q- 

The value of (yFeji) in Equation (7) is within the range of possible SNII iron yields calculated by 
Woosley & Weaver (1995); this provides some confidence in the veracity of other parameters determined 
from the equations above. However, the specific supernova rate r\n — 0.026 is almost four times the value 
expected for a Salpeter IMF, r\ st d = 0.00681 (x = 1.35, m u = 100 M and m e = 0.08 M ). Values of rj u 
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computed with single power law IMFs are sensitive to both the slope x of the IMF and the mass limits m u 
and me. The variation of 7777 and (3 with x, m u and me for several power law IMFs is illustrated in Figure 2. 

It is difficult to assign errors to the quantities evaluated in Equations (7) - (9) because of the many 
uncertainties involved and the possibility of systematic errors. If all the quantities in Equation (8) are 
skewed to their limits of uncertainty in an effort to minimize 77//, it is possible to achieve a close agreement 
with a Salpeter IMF. For example, if T S i, g = 0.0225, R = 0.92, M/L B = 9, (y S i,n) = 0.32 (model WWlex 
from Loewenstein & Mushotzky 1996), and (3 = 0.3 then rju = 0.0068, close to the Salpeter value. Of 
course it is most unlikely that Nature would conspire in this manner. 

Also shown in Figure 1 are the fraction of all iron created by SNIa: 

p _ niaUFeJa 

Ia n Ia y Fe j a + (1 - (3)- 1 {M/L B )r] II {y FetII ) 

and the fraction of all iron in the hot ICM gas that originated in SNIa: 

p n Ia y Fe j a (l - Gt) 

Ia ' 9 ~ n Ia y Fe . Ia (l - G„) + (1 - (3)^{M/L B ) m i{y Fe , n ){l - F#) ' 

Provided SNu(t„) is limited to values of interest, log[SNu(t„)] < — 1, Type Ia supernovae can contribute to 
the stellar enrichment without substantially changing r/n or (y Fe ,n)- 

Renzini et al. (1993), Ishimaru & Arimoto (1997), Renzini (1997) and Wyse (1997) have discussed 
the X-ray data for rich clusters in detail. These authors prefer models in which the overabundance of 
iron relative to that expected with a normal Salpeter IMF is due to a large additional iron contribution 
from Type Ia supernovae, corresponding to Fj a ~ 0.5 in Figure 1. Such a model would be very similar to 
the enrichment history of our own Galaxy. In particular, Renzini et al. (1993) suggest a model for rich 
cluster enrichment in which ~ | of the iron is produced by SNIa with the remaining ~ | coming from 
SNII. If only half of the iron in rich clusters has originated in SNIa, then the number of SNIa per Lb is 
nia = A/jo/ £b ~ 0.025. But such a large number of SNIa is incompatible with iron and silicon abundances 
in our calculated model of a single large elliptical, discussed below, where we find m a = Nia/C-B ^ 0.001. 
For the simple power law model for SNu(t) we adopt in Figure 1 (p = 1, t* = 2 Gyrs, t n = 13 Gyrs), the 
current SNIa rate would need to be very large, SNu(i„) > 0.6, for Type Ia supernova to supply half of the 
iron. This is about ten times greater than the SNIa rate estimated from observations of bright ellipticals. 

An alternative explanation of the iron excess in rich clusters relative that produced by a Salpeter 
IMF (x, m u , me — 1.35, 100, 0.08) is to invoke a flatter IMF having a larger fraction of massive stars 
and associated SNII. This is the interpretation preferred by David (1997) and Gibson, Loewenstein & 
Mushotzky (1997). For consistency both (3 and r/ii must be adjusted when the IMF slope is changed; the 
value of the mass return fraction (3 = 0.3 used in Equation (8) is based on the Salpeter IMF. Logically 
compatible parameters can be found by simultaneously solving Equation (8) for nn{(3) with the parametric 
variations rju (x) and (3(x) plotted in Figure 2. The result of this joint solution is that Equation (8) now 
gives nii = 0.016 with j3 = 0.55 for (x, m u , m t = 1.00, 100, 0.08) or n n = 0.018 with (3 = 0.51 for (x, m u , 
m e = 0.94, 40, 0.08). (Such large values of the mass return (3 from evolving stars may be unrealistic since 
several generations of stars must have formed near time t* , each consuming some of the metal-enriched 
gas ejected from previous generations of stars.) With these flat IMFs almost all of the iron in rich clusters 
originates in SNII and the low SNIa rates are consistent with those observed in massive ellipticals. IMFs 
flatter than Salpeter are also supported by the large mass of oxygen in rich clusters; Gibson, Loewenstein 
& Mushotzky (1997) claim that the observed O/Si ratio can be produced by SNII alone with little or no 
contribution from SNIa. 
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When silicon and iron abundances in the ICM are considered together with abundance constraints set 
by individual massive ellipticals, as we have done here, the metal enrichment of rich clusters by early type 
galaxies is shown to be almost exclusively due to Type II supernovae with a negligible contribution from 
Type la supernovae. Evidently, star formation conditions in rich clusters are quite unlike those that prevail 
in our own Galaxy. 

The controversy concerning the relative enrichment from Type II and Type la supernovae in rich 
clusters is further confounded by the gas dynamical solutions we discuss below. We describe a variety of 
supernova enrichment histories for massive elliptical galaxies in which the hot interstellar gas is constrained 
to evolve toward the radial distribution of gas density, temperature and metal abundance observed today. 
In particular we show that the production of SNII using a Salpeter IMF is sufficient to achieve simultaneous 
agreement with n(r), T{r) and zp e {r) indicated by X-ray observations of massive ellipticals. Gas dynamic 
solutions based on IMFs flatter than Salpeter are generally incompatible with the observed abundance 
variation z Fe (r). The SNIa rate must also be low. If more supernovae of cither type are involved, as 
indicated by the large r\u for rich clusters (Equation 8), the iron and silicon abundances in the models 
greatly exceed those observed in massive elliptical galaxies. 

3. IRON AND SILICON IN NGC 4472 

3.1. Gas Dynamical Models for NGC 4472 

In a recent paper (Brighenti & Mathews 1998b) we discuss in detail the basic assumptions and 
equations used in our gas dynamical models that simulate the evolution of hot interstellar gas in elliptical 
galaxies. Since the models we discuss here are very similar to those in that paper, we provide only a brief 
summary. 

Our ID spherical calculations begin with an overdensity perturbation in a flat cosmology having an 
amplitude designed to produce a mass concentration similar to that of NGC 4472 after a few gigayears. The 
flow of dark and baryonic matter far from the center of the perturbation follows the self-similar solution 
described by Bertschinger (1985). An outward moving turn-around radius defines the instantaneous locus 
where the cosmic flow velocity vanishes. The "secondary infall" occurs within the turn-around radius, 
allowing baryonic and dark matter to collect near the origin. The central accumulation of dark matter grows 
from the inside out; although the collisionless dark fluid interpenetrates, after passing through the origin 
it continues to oscillate, spending most of its time at large radii. The net effect is that a quasi-stationary 
core of dark matter grows outward from the origin and is matched to the outer Bertschinger inflow in a 
manner that preserves the total mass of dark matter. However, we replace the inner power law core of the 
self-similar Bertschinger solution with a (less peaked) dark halo having an NFW profile as determined by 
Navarro et al. (1996) with full three-dimensional N-body calculations. The dark matter mass of NGC 4472 
inferred by X-ray data can be fit reasonably well with an NFW dark halo of mass Mh = 4 x 10 13 M Q within 
the current virial radius assuming D = 17 Mpc for the distance to NGC 4472. The cold baryonic gas also 
participates in the Bertschinger flow, but deep within the turnaround radius it passes through an accretion 
shock, compresses and is heated to approximately the virial temperature of the galactic dark halo. At time 
t* when enough baryonic matter has collected within the shock, some of which has radiatively cooled, we 
form the de Vaucoulcurs stellar configuration having a total mass = 7.26 x 10 11 M Q appropriate to 
NGC 4472 and corresponding to a stellar mass to light ratio M*/Lb = 9.20 (van der Marel 1991). The 
stellar galaxy is constructed by removing a baryonic mass equal to M* from gas within the accretion shock 
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r s h(t*) to form the stars. The gas density in the remaining gas is reduced in proportion to its density just 
before t*. 

We assume the SNII energy is released immediately at time t* . After gas has been removed to form the 
stars, the SNII energy is assumed to be evenly deposited (masswise) to remaining gas within the accretion 
radius or within some other specified radius. Removing gas interior to r s h{t*) and heating the remaining 
gas introduces a transient in the flow, but after <~ 1 Gyr gas moving within r s h becomes subsonic and 
approaches hydrostatic equilibrium. The subsequent time dependent gas dynamics within and around 
the galaxy are followed in full detail. The potential of the dark matter continues to evolve according to 
the NFW-Bertschinger prescription and the baryonic shock grows in radius. The flow of these two fluids 
is solved simultaneously using ID Eulerian hydrodynamics on a logarithmic grid. The equations that 
describe this flow are described in detail in our recent papers (e.g. Brighenti & Mathews 1998b). Inside 
the de Vaucouleurs core new gas is provided by mass loss from stars as they evolve off the main sequence 
according to some assumed stellar IMF. Most of the gas within the accretion shock radiates thermal X-rays 
corresponding to a temperature T ~ 10 7 K similar to the virial temperature in the total potential. After 
time t* additional energy is supplied to the gas by the dissipation of the orbital energy of mass-losing 
stars and by Type la supernovae. The stellar velocity field is found from solutions of the Jeans equation 
in the combined stellar and dark matter potential, although the evolved solutions are (fortunately) rather 
insensitive to the precise variation of the stellar velocity ellipsoid with galactic radius, as characterized by 
the stellar temperature T*(r). This galaxy formation model is obviously oversimplified, but it properly 
conserves both dark and baryonic matter and provides a realistic template for our gas dynamical models. 

In a typical calculation the shock radius at time t* is r s h{t„) <~ 200 kpc, considerably larger than 
the present size of the stellar system for NGC 4472 (effective radius r e — 8.57 kpc, maximum stellar 
radius r*i = 100 kpc). It may seem inconsistent to introduce SNII heating and metal enrichment within 
a radius r s h{t*) that is much larger than the currently observed size of the stellar galaxy. This apparent 
inconsistency can be justified because star formation must have occurred in a somewhat larger region than 
is occupied by old stars in NGC 4472 today; in order to form the de Vaucouleurs profile in NGC 4472 
by violent relaxation, it is necessary that the pre-galactic objects be largely stellar before the galaxy is 
assembled. 

The total energy released by SNII explosions at time t* is En = e sn r]jjM^E sn where E sn = 10 51 ergs, 
■A/jj = r/u-M* is the total number of SNII produced and e sn is the efficiency that SNII energy is converted 
to thermal energy in the ambient gas, the rest is lost to radiation. For a reference value of r\n we use the 
value computed with a Salpeter IMF (x, m u , mi = 1.35, 100, 0.08), rj st d — 6.81 x 10~ 3 . For simplicity, 
the stellar mass in these relations (and in establishing the fixed stellar potential) is set to the current 
value M* = M*(t n ) rather than its value at t*, M*/(l — [3). In view of the uncertainties involved in all 
these parameters, in the gas dynamical models we evaluate En = r/j/M* (f„)10 51 as if e sn = 1 and vary 
r/jj, seeking results that best match current observations. Each Ejj determined in this way corresponds 
to a range of M*(t*)E sn for each assumed gas heating efficiency e sn — Eji /r)nM*(t*)E sn < 1. Metal 
enrichment by SNII accompanies the energy deposition within the accretion shock radius at t* . Altogether 
(yFe.ii)vn^* solar masses of iron are available for enriching stars and gas within r s h(t*). If zp e ,* is the 
assumed stellar iron abundance, a mass ZF e ,*M*/lA of iron produced by SNII is allocated to the newly 
formed stellar system at time t* and the remainder is introduced into the gas within the shock radius in 
proportion to the local gas density. 

In summary, the important galactic parameters that determine the early star formation and SNII 
enrichment in rich clusters and individual massive ellipticals are (i) the total energy released in SNII 
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explosions En , (ii) the number of SNII per solar mass ?]n , and (iii) the time of star formation t* . 

3.2. The Standard Model 

In our standard or reference model for the evolution of hot interstellar gas in NGC 4472, the galaxy 
forms in a flat universe (CI = 1) with global baryon fraction Ctb = 0.05. The specific SNII rate is 
f]n = Vstd = 6.81 x 10~ 3 . The specific stellar mass loss rate for the Salpeter IMF a*(t) — dlogM*/dt is 
well fit with a power law a*(t) = 4.7 x 10 _20 (i/i s )~ 126 s _1 where t n = 13 Gyrs and and t s = t n — i* is the 
current age of the stars. Both stars and SNII are assumed to form at time t* = 2 Gyrs. The total number 
of SNIa, A/j , is an integral over SNu(i) which is assumed to vary as a simple power law parameterized with 
SNu(i„) = 0.015 SNu and p = 1. Type la supernovae are assumed to begin at i* and continue until t n = 13 
Gyrs. Except for the slightly lower SNu(t„), all parameters are identical to the standard reference model 
discussed in Brighenti & Mathews (1998b). In addition, we assume for the standard model that no iron or 
silicon produced by SNIa is used in stellar enrichment, so the stars have abundances proportional to SNII 
yields. 

The first requirement for an acceptable model is that the chosen parameters, when used in a gas 
dynamical calculation, adequately reproduce the radial variation of interstellar gas density and temperature 
observed in NGC 4472 today. In Figure 3 we plot with solid lines the density and temperature profiles 
in the hot interstellar gas for the standard model at time t n = 13 Gyrs. The overall agreement with the 
observed gas temperature and density is satisfactory but not perfect. The excess gas density in the model in 
r < r e = 8.57 kpc is a classical artifact of galactic cooling flow models; we believe that it can be mitigated 
by galactic rotation and low mass star formation but this has not yet been adequately demonstrated (see 
Mathews & Brighenti 1998c for a brief review). At larger radii, r > 10 kpc, the computed density is 
slightly low. This latter discrepancy may be due to the disturbed nature of NGC 4472 at large radii where 
it appears to be interacting with ambient gas in the Virgo cluster (Irwin & Sarazin 1996) although the 
azimuth-averaged outer gas density profile in NGC 4472 is typical of the X-ray structure in other large 
ellipticals (Mathews & Brighenti 1998b). In any case, if En is slightly lower than the Salpeter value, the 
standard model can be adjusted to fit the observations almost perfectly at large r (Brighenti & Mathews 
1998a). 

For the models discussed here we assume a simple flat universe, CI = 1, but it is also possible to find 
satisfactory and very similar dynamical models for NGC 4472 using other cosmologies. In the Appendix we 
briefly describe gas dynamical models appropriate for an open universe (Cl = 0.3) and a low density flat 
universe (SI — 0.3, Cl\ = 0.7). 

A fully successful model must also match the metal abundances and abundance gradients currently 
observed in NGC 4472. Unfortunately, supernova yields and observed abundances are both uncertain so 
there is some flexibility in making this fit. The iron abundance in NGC 4472 has been observed extensively 
with ROSAT and ASCA by many observers using many different data reduction procedures. Observed iron 
abundances expressed in solar units can be a source of some confusion if authors do not explicitly note which 
iron abundance they have assumed for the sun, photosphcric or meteoritic. If the lower meteoritic value is 
more correct, as we assume here, iron abundances relative to the solar photospheric value are too low by 
1.44 (Anders & Grevesse 1989). In the following discussion, we consider only observational sources for which 
the absolute solar abundance is specified and, if necessary, we convert observed iron abundances cited in 
the literature to meteoritic solar. Unfortunately, observational determinations of the iron abundance in the 
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hot gas also seems to depend on the data reduction procedure used. Some of the difficulty in deriving hot 
gas abundances arises due to the presence of additional X-ray emission unrelated to the hot interstellar gas. 
Harder radiation (E > 2 keV) thought to be stellar in origin must be allowed for in fitting the spectrum 
and deriving abundances. Observational determinations of the iron abundance appear to depend on the 
procedure used in allowing for the hard radiation. For example, Buote and Fabian (1998) find that the iron 
abundance in NGC 4472 determined with a single temperature plasma increases by a factor of 4.7 when a 
two-temperature fit is used instead. Possible inadequacies in treating the Fe L transitions in the adopted 
plasma code may also influence the resulting value for the observed iron abundance. 

Seeking some consistency, we consider here several published globally averaged iron abundances 
determined with detectors on the ASCA satellite. Using the Raymond-Smith plasma code with a hard 
component in the spectrum, Arimoto et al. (1997) find zp e = 0.47 ± 0.09 solar in the hot interstellar gas 
in NGC 4472 within 3 arc minutes (15 kpc) of the galaxy center. Matusmoto et al. (1997) use a fixed 
hard spectrum typical of low mass binary stars and find zp e = 0.56 ± 0.09 solar for NGC 4472. Finally, 
Matsushita (1997) gives ZFe = 0.63 ± 0.08 solar for the iron abundance in NGC 4472. All these values 
are consistent at about half solar, but Buote & Fabian (1998), also using ASCA data with the MEKAL 
code, derive two quite different abundances: zp e = 0.37 solar (single temperature model: T = 0.97 keV) 
and ZFe = 1.70 ± 0.46 solar (two temperature model: T\ = 0.76 keV, T 2 = 1.48 keV). All these values are 
emission-weighted (oc p 2 ) and therefore apply mostly to denser interstellar gas located closer to the galactic 
center. In addition to these globally averaged values, Matsushita (1997) finds a negative iron abundance 
gradient in several bright ellipticals including NGC 4472; this gradient will play an important role in fitting 
to abundances computed in our gas dynamical models. Because of the additional spatial information 
provided, we adopt Matsushita's iron abundances for comparison with our models. Finally, Matsushita 
(1997) and Arimoto et al. (1997) find a global silicon abundance z$i ~ 0.5 solar for NGC 4472. 

Figure 4 illustrates the distribution of iron and silicon in the hot gas surrounding NGC 4472 after 
evolving to the present time t n — 13 Gyrs. This calculation is based on supernova parameters for the 
standard model - SNu(i„) = 0.015, p = 1 and mi = 6.81 x 10 -3 - with supernova yields: (yFe.ii) — 0.14 
M , (ysiji) = 0.133 M , y Fe ja = 0.744 M Q , and y Si ,i a = 0.158 M . Following Arimoto et al. (1997), we 
assume a power law gradient for the stellar metal abundance 

z Fe ,* = 0.675(r/r e )-°- 207 

which is based on the observed spatial variation of the Mg2 photometric index. With this variation the 
mean stellar Fe/Mg ratio is half solar, (zF e ,*) / (zMg,*) — 0.5, (in solar units) as suggested by Trager (1997) 
where (zmj,*) = 1-385 solar (Arimoto et al. 1997). Since SNIa do not enrich the stars in the standard 
solution, the stellar silicon abundance is proportional to SNII yields, z$i,* — ZFe,*(ysi,n)/ (yFe,ii)- 

In the standard model all SNII occur as the stellar system is created at time t* = 2 Gyrs and the stars 
are enriched with enough iron and silicon from SNII to match the mean stellar abundances (zpe.*) and 
(zsi,*)- All additional iron and silicon from SNII is distributed at time t* to gas within the accretion shock 
radius r s h(t*) = 200 kpc in proportion to the local gas density. The energy released by SNII is sufficient to 
temporarily reverse the secondary inflow, causing the contact discontinuity that defines the enriched region 
within r s h(t*) to expand to 400 kpc at the present time as shown in Figure 4. The uniformity of zp e and 
zsi in the "plateau" region visible in Figure 4 between the outermost stars and the contact discontinuity, 
100 - 400 kpc, is an artifact of the uniform deposition of SNII metals at time t*. Most of the iron mass 
is in this plateau region as shown in the plot of Mp e {r) in Figure 5. At time t n the interstellar gas is 
slowly flowing inward within the galaxy r ^ 100 kpc, where the early SNII enrichment is diluted by new 
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gas contributed by stellar mass loss. But this dilution is more than compensated by additional iron and 
silicon contributed by SNIa and stellar mass loss within the galaxy. The total iron abundance distribution 
from all of these sources, shown with a solid line in Figure 4, is seen to agree reasonably well with the 
observations of Matsushita (1997). Emission-weighted abundances based on this standard model depend 
on the radius considered. For the entire region within the current accretion shock radius r s h(t n ) = 1000 
kpc the mean abundances in the gas (in solar units) are (zp e ) — 0.45 and (zsi) — 0.52; within the current 
contact discontinuity at 400 kpc the mean gas abundances are (zp e ) = 0.67 and (zsi) = 0.77; within 150 
kpc the abundances are even higher, (zFe) = 0.76 and (zsi) = 0.86. These latter values are slightly higher 
than the global iron and silicon abundances reported by Matsushita (1997); this excess probably arises 
because the gas density in our models is slightly larger and more centrally peaked than that observed in 
NGC 4472 (Figure 3). The approximate constancy of (zsi/ 'zsi,®) I '( z Fe/ ' z Fe,o) ~ 0.6 with galactic radius 
and its value intermediate between SNII and SNIa yields ({ysiji) / il/Feji) ~ 0.8 and ysi.ia/lJFeja ~ 0.2 
respectively) are both attributes of recent X-ray observations of NGC 5846 by Finoguenov et al (1998). 

For any given gas dynamical solution described by n(r) and T(r) at time t n , an infinite family of 
interstellar abundance distributions can be generated by changing the specific number of SNII rju (or the 
supernova yields) while keeping the total SNII energy En = e sn 7]uM*E sn unchanged. The constancy of 
En is ensured by varying the gas heating efficiency factors e s „ or E sn to compensate for changes in rju. 
When following this procedure, we were surprised to discover that the iron abundance distribution shown in 
Figure 4 disagrees significantly with Matsushita's observations when rju is only varied by ~ 10 percent from 
that used in the standard solution. For example, if r)n — 1.25r] st d the iron abundance in the plateau region 
(100 - 400 kpc) increases to zpe/zFe.Q — 0-72 and the central value is nearly 1.1. These iron abundances 
are significantly higher than Matsushita's observed values, particularly in the outer galaxy r > 20 kpc. The 
sensitivity of zp e (r,t n ) to r\n can be understood because of the large amount of iron and silicon required 
to enrich the galactic stars. The global iron abundance in the gas just after SNII enrichment, identical to 
that in the plateau region, is zp e = (Mpe.ii — Mp e ^)/M gas where Mpe,n is the total iron produced by 
SNII. Since the large total mass of iron in stars Mp e ,* is held fixed as r\n is varied to create a family of 
dynamically identical enrichment models, small changes in i]n and therefore Mpeji correspond to rather 
large changes if zp e since Mp e ,ii is not much larger than Mp e ,*- ln this sense our models require a high 
degree of regularity in rju in order to match the negative gaseous iron abundance gradients typically 
observed in giant ellipticals. This sensitivity to rju may be an indication that our model is too simple or 
that some additional enrichment process has been overlooked. 

Although our standard model, as illustrated in Figures 3 and 4, is successful in approximately 
reproducing all relevant observations of NGC 4472, the value of rju that we have used corresponds to a 
normal Salpeter IMF and is therefore several times less than values required to enrich massive clusters of 
galaxies (Equation 8). In spite of this important difference in rju, the gaseous iron and silicon abundances 
in our standard model within the current accretion shock radius r s h(t n ) = 1000 kpc, (zp e ) = 0.45 zp e , Q and 
{zsi) = 0.52 ZFe,®, are both rather similar to typical abundances in the ICM of rich clusters, ~ 0.4 solar. 
The ratio of total iron mass in both stars and gas to the total baryonic mass within r s h(t n ) corresponds to 
an iron abundance of ~ 0.4 in solar units, again very similar to cluster values. This apparent agreement 
is likely to be just a coincidence, however. The important distinction between our model for NGC 4472 
based on a Salpeter rju and higher values of rjn required to understand cluster ICM abundances is in 
the much lower total gas mass fraction currently present in galaxies as compared to clusters. Within any 
radius of interest out to the current accretion shock radius, the ratio of gas mass to total baryonic mass 
in our standard model is much lower than that found in rich clusters. If the ICM in rich clusters contains 
additional, unenriched primeval gas at larger radii, the enrichment of cluster gas by ellipticals similar to our 



- 15 - 



model for NGC 4472 is even more dramatically inadequate. Therefore, our models support the contention 
made by Renzini (1997) and David (1997) that it is impossible to build presently observed rich clusters by 
merging hot gas in group-dominant massive ellipticals having properties similar to our standard solution for 
NGC 4472. 

We have made additional calculations using standard model parameters but assuming that the SNII 
energy and enrichment occur within radii that are not equal to r s h(t*). When the enrichment region is 
smaller than r s h (t* ) , the plateau region is smaller but its iron abundance is increased so that it is no longer 
possible to match the current negative iron abundance gradients observed in NGC 4472 and other similar 
massive ellipticals. To correct for this excess iron, the value of r\n would need to be less than r\ s td-, diverging 
further from typical cluster values. If the size of the enrichment radius is much larger than r s h{t*), then the 
dynamical time for other stars spatially associated with the enriching SNII to arrive at the outer radius of 
the currently observed stellar system in NGC 4472 exceeds the time available t n — i*. For example, stars 
that form at t* beyond r w 700 kpc in the standard solution can never enter the stellar part of the galaxy 
(r < 100 kpc) by the present time t n = 13 Gyrs; any SNII enrichment beyond this radius is unrelated to 
stars within NGC 4472. In any case, even if we use a SNII energy and enrichment radius equal to the 
maximum allowed size at i*, <~ 700 kpc, the current abundances of iron and silicon are still too large if the 
cluster value of r\n is used. 

3.3. The Cluster Model 

To emphasize the disparity between optimum values for the specific number of SNII in rich clusters 
and in single bright ellipticals, we briefly discuss a "cluster" model for NGC 4472 in which r\n has a larger 
value appropriate for rich clusters. Following our previous discussion in §2, we assume rju = 2.i&rj st d 
with (3 = 0.55 and a slightly larger mean SNII iron yield, (yFeji) — 0.158. All remaining variables and 
the gas dynamical solution are identical to those in the standard model previously discussed. Larger rju 
corresponds to a lower heating efficiency e sn , keeping Ejj and the associated gas dynamics unchanged from 
the standard solution. 

After evolving to the current time, the iron and silicon distributions in this model, shown in Figure 6, 
are clearly unable to account for gas abundances observed in NGC 4472. Iron and silicon within the galaxy 
are dominated by the inflow of SNII-cnriched gas from outside the stellar system. The total iron abundance 
within the optical galaxy (r < 100 kpc) is much higher than observed abundances and its radial gradient 
is positive throughout the galaxy. This model is an example of the sensitivity of ZFe(r) to the parameter 
rjji. Increasing rjjj by a factor of 2.36 causes the gas abundance in the plateau region to increase to 3.35, 
almost 22 times greater than the abundance in this region in the standard solution. As before most of the 
iron mass in the gas is contained in this distant plateau region (see Figure 5). 

We explored several variants of the "cluster" model in an attempt to reconcile cluster values of r]u 
and the enrichment history of large ellipticals. The high abundances in the "plateau" region just beyond 
100 kpc can be reduced if the radius of SNII enrichment at time i* is increased or if this gas is (rather 
arbitrarily) mixed and diluted with primordial gas at larger radii. However, as discussed above, such models 
are unrealistic since the dynamical times at ~ 700 kpc for either of these variant models is too long. Low 
mass stars accompanying metal-producing SNII at these radii will not have joined the galaxy by the present 
time. We have also sought solutions in which the galaxy formation time t* is 3 or 4 Gyrs in an attempt 
to lower the abundances by mixing SNII ejecta with larger masses of unprocessed gas inside the accretion 
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shock at these later times, but unrealistically high iron abundances are still present. 

We are faced with the unexpected result that the amount of metals required to fit abundances in rich 
clusters are incompatible with those observed in individual ellipticals. This is very curious since bulge 
dominated E and SO galaxies are generally thought to be the primary source of metals in the ICM of rich 
clusters. 



3.4. The la Model 

We now explore the possibility that SNIa are important contributors to the iron abundance in rich 
clusters and that bulge dominated galaxies are the primary sites of these SNIa explosions. This can be 
accomplished with the "la" model for NGC4472 for which nj a = 0.011, corresponding to Fi a = 0.22 and 
Fia.g — 0.25 in Figure 1. To achieve this high SNIa rate we choose SNu(t n ) = 0.1 and p — 1.5; this is the 
same slope as that of Ciotti et al. (1991), but our SNu(t n ) is only 70% of their value. All other parameters 
are identical to those in the standard model. 

Since the energy deposited by SNIa in this model exceeds that of the standard model, the hydrodynamic 
solution is also different. The gas density and temperature for this model after time t n = 13 Gyrs are 
shown with dashed lines in Figure 3. The agreement with the observed gas density in NGC 4472 is good, 
but the gas temperature has a minimum near 60 kpc that is not observed. This curious feature is a 
long-lived relic resulting from gas flows driven by SNIa and SNII heating at times just after time t*. In 
earlier theoretical models for the evolution of hot gas in ellipticals in which stars are assumed to be the only 
source of interstellar gas, such a large energy deposition by SNIa would have driven a strong galactic wind. 
For the model we discuss here, where the young galaxy is surrounded by cosmic gas converging toward the 
overdensity perturbation (secondary infall), early SNIa-driven winds are suppressed by the inertia of this 
ambient gas and only a modest outward redistribution of the gas occurs. 

Apart from the small temperature minimum in Figure 3, the general trend in the gas temperature for 
this model is a reasonably good fit to gas temperatures observed in NGC 4472. However, the current iron 
and silicon abundances in the hot gas for this model, illustrated in Figure 7, are clearly at variance with 
X-ray observations. In this SNIa-dominant model the interstellar enrichment of iron and silicon is caused 
almost exclusively by SNIa, raising the iron abundance to > 6 times solar throughout the optical galaxy. 
The shallow positive gradient dzFe/dr > for r < 100 kpc is due to a small residual gaseous outflow 
caused by the large energy released by SNIa inside the galaxy. The global value of the iron abundance 
and its radial gradient are both too large. This model is rather similar to the earlier result of Loewenstein 
& Mathews (1991) in their study of simpler galactic cooling flows without cosmological secondary infall; 
they also found iron abundances far in excess of those observed when SNIa are important contributors 
to the overall galactic enrichment. It is also significant that the silicon abundance in this model is much 
higher than the global value observed in the hot interstellar gas of NGC 4472, (zsi) — 0.5 solar. Most of 
the excess iron and silicon inside the galaxy in this model comes from Type la supernovae with a smaller 
contribution from stellar mass loss. Because of the low-level, subsonic outflow driven by SNIa, iron and 
silicon introduced into the gas by SNII at £* are unable to enter the interstellar gas interior to ~ 100 kpc. 

We have also explored additional high nj a models in which the SNIa varies like a step function rather 
than a power law. By concentrating most of the SNIa energy at early times, we hoped that the large amount 
of metals produced then would be expelled in a powerful galactic wind. If so, the low gas abundances 
observed today within ~ 100 kpc might be made consistent with large SNIa enrichment at early times. For 
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example, we considered a model with n/ a = 0.05, with most of the SNIa exploding at very early times: 
SNu = 24.86 for 2 < t < 4 Gyr, and SNu = 0.03 for t > 4 Gyr. However, with this sort of model the gas 
dynamics were altered so that we were unable to match n(r) and T(r) currently observed in NGC 4472. 
The iron and silicon abundance distribution are also strange in these models at the present time, with an 
enormous metal enrichment occurring just outside of the stellar part of the galaxy. 

4. Conclusions and Final Remarks 

When we began this study of the enrichment history of massive ellipticals and rich clusters, we quite 
mistakenly anticipated that it would be possible to find a common set of parameters that would account for 
evolutionary enrichment on both galactic and cluster scales. Gas phase iron abundances in bright ellipticals 
range from 0.1 — 1 solar, with NGC 4472 lying near the upper limit of this range, but the average may 
not differ greatly from the mean global iron abundance in rich clusters, ~ 0.4. Rich clusters typically have 
higher gas fractions than massive ellipticals, but this might result in part from supernova-driven outflows 
from the shallower potentials of member ellipticals. Unlike the observational information available for hot 
gas in group-dominant ellipticals, our computational models are closed boxes since we can computationally 
follow gas-phase supernova enrichment products into regions of very low gas density and emissivity where 
observations are difficult or impossible. Perhaps, we supposed, the large number of past supernovae required 
to produce metal abundances in rich clusters were at one time also present in giant ellipticals but some 
of these metals was dispersed by galactic winds into the low density environment, escaping observational 
detection. If all of this were true, it might be possible to merge small elliptical-dominated groups together 
with their wind-enriched gaseous environments and construct rich clusters. In developing this hypothesis, 
we adopted the commonly held assumption that virtually all metal enrichment in rich clusters is due to the 
stellar spheroidal component, i.e. in E and SO member galaxies. 

As discussed in our recent paper (Brighenti & Mathews 1998b), the observed hot gas density and 
temperature distributions in any massive elliptical can be fit with a range of gas dynamical models having 
different numbers of Type II supernovae Nji, represented here by the SNII number per unit stellar mass, 
rjij = Njj/M*, or the total energy released by all SNII, En = e sn rjuM*E sn . For a given set of cosmological 
parameters, hydrodynamical models for the evolution of hot interstellar gas are degenerate in the sense that 
a range of values of r\n (associated with different IMFs) and therefore Ejj can produce similar interstellar 
density and temperature distributions at the present time provided the time of galaxy formation £* is also 
adjusted. If the galaxy is assumed to form at a later time, more gas is available near the central overdense 
region and a larger value of Ejj generates approximately the same specific thermal energy; the subsequent 
evolution of the gas arrives at similar n(r) and T(r) distributions at the present time. This approximate 
degeneracy depends only weakly on our assumption that the SNII energy is deposited within the accretion 
shock radius at time t*, r s h{t*)- As we have discussed already, the extent of the SNII enrichment region rjj 
in massive ellipticals cannot be greatly different from r s h{t*)- If rjj < r s h{t*), then even lower values of X]n 
arc required to reproduce abundances observed today, increasing the disagreement with the higher cluster 
values for r\n. If rjj > r s h(t*), then the dynamical time also increases and low mass stars associated with 
the SNII enrichment would not have enough time to enter the massive optical galaxy observed today. 

The important parameters that describe our evolutionary models are (i) the total energy released by 
Type II supernovae, Ejj = e sn r]iiM*E sn which must be adjusted to achieve agreement with n(r) and 
T(r) in the hot interstellar gas, (ii) t* the time of star formation, (iii) rju which controls most of the gas 
enrichment, and (iv) the SNIa rate SNu(t) that supplies much of the interstellar iron in the central galaxy, 
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r < 50 kpc. In seeking satisfactory gas dynamical models with no attention to abundances, the main 
parameters are En and i*. We have found that models for the large elliptical NGC 4472 agree with n(r) 
and T(r) and the abundance distribution zp e (r) if stars form at = 2 Gyrs with a Salpeter IMF for which 
Vn = Vstd = 6.81 x 10~ 3 assuming a heating efficiency e sn = 1. If we wish to increase r) U and therefore 
En toward the global values indicated by rich cluster abundances, ~ 2 — 3n s td, t* must be increased 
significantly. However, such an increased i* is not possible since observations indicate that many luminous 
elliptical galaxies are very old: the color-magnitude diagram (Bower, Lucey, & Ellis 1992), the small scatter 
in the Mg2 — a relation (Bender, Burstein, & Faber 1993), passive evolution of the fundamental plane (van 
Dokkum & Franx 1996), etc. With i* constrained to values 2 Gyrs, and En constrained for a match 
currently observed n(r) and T(r) distributions, rjn can be increased toward the cluster value (~ 2.6rj st d) 
only by decreasing the gas heating efficiency e sn with En and held fixed. This is the procedure we have 
considered with the "cluster" or high-77/7 model for NGC 4472. However, as rjii is increased above the 
standard Salpeter value, the metal enrichment of the interstellar gas at the present time quickly rises far 
above abundances observed in NGC 4472. The global iron and silicon abundances and the negative radial 
iron gradient dzp e /dr in NGC 4472 (and other large ellipticals) can be matched only with a very narrow 
range of rjn, all lower than cluster values. Similar restrictions also apply to the total number of Type la 
supernovae which we assume are also more frequent in the early universe. 

We arrive at the following conclusions: 

(1) Type II supernova production with a Salpeter IMF is sufficient to explain all important hot gas 
observations in bright ellipticals: the radial distribution of density, temperature and iron abundance as well 
as the global iron and silicon abundances. A small additional contribution of iron and silicon from Type la 
supernovae, consistent with currently observed SNIa rates in elliptical galaxies, can also be present. Unlike 
the enrichment history of our Galaxy, where iron production by Type II and Type la supernovae have been 
comparable, in ellipticals the total iron and silicon enrichment by Type la supernovae must be much less 
than that from SNII. 

(2) The gas dynamic and enrichment history of massive ellipticals establishes upper limits on the total 
number of both Type II and Type la supernovae events. Large numbers of Type II supernovae associated 
for example with flat IMFs, often invoked to account for the iron observed in rich clusters, produce too 
much gas-phase iron and silicon in models of single massive ellipticals. If the total number of Type la 
supernovae is large enough to dynamically influence the interstellar gas in these galaxy models, the current 
metal enrichment greatly exceeds the gas phase iron and silicon abundances observed. 

(3) For fixed total SNII energy En, the iron and silicon enrichment is very sensitive to the specific 
number of Type II supernovae r\n involved. This sensitivity may indicate that our galaxy enrichment model 
is oversimplified or it may be related to the wide spread in gas phase abundances observed in ellipticals and 
group-dominant ellipticals (see Fig. 6 of Renzini 1997). 

(4) We agree with David (1997) and Renzini (1997) that it is not possible to form presently 
observed rich clusters by simply combining stars and hot gas from a large number of presently observed 
elliptical-dominated groups. 

(5) In view of our model calculations for the evolution and enrichment of interstellar gas in massive 
ellipticals, neither of the currently-discussed explanations of cluster metal abundances is particularly 
attractive: (i) IMFs flatter than Salpeter, and (ii) substantial iron enrichment from Type la supernovae, as 
in our Galaxy. 
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These conclusions could be changed if some of the theoretical and observational parameters we have 
used are in error by large amounts - for example if observed metal abundances in the hot interstellar gas of 
ellipticals are underestimated by factors of 5 or more - but this seems unlikely. 

Alternatively, the different enrichment histories of large ellipticals and rich clusters could be reconciled 
if the cluster gas is enriched by another source not involving massive elliptical galaxies. Arnaud et al. 
(1992) show that the total mass of hot gas in clusters correlates nicely with the total optical luminosity of 
E and SO galaxies within 3 Mpc of the cluster centers, but there is no similar correlation with the total 
luminosity of the spiral galaxy component. Since global gas phase iron abundances are fairly constant 
among clusters, a correlation with gas mass is also a correlation with the total mass of iron, most of which 
is in the gas phase. Arnaud ct al. (1992) therefore concluded that E and SO galaxies are the only significant 
sources of metal enrichment in the hot cluster gas. 

However, cluster enrichment may be more complicated than this. For example, the role of disk galaxies 
in ICM enrichment may be more important than is currently thought. Most low redshift rich clusters 
contain about twice as many SO galaxies as ellipticals (Dressier 1980). The fractions of SO and spiral 
galaxies in clusters are inversely related while the ratio of the number of SO plus spiral galaxies to the 
number of ellipticals is rather constant. This suggests that SO galaxies may be descended from spirals. 
This is supported by the observation that richer clusters tend to have higher SO/spiral ratios. Even more 
relevant, cluster observations at z ~ 0.5 reported by Dressier et al. (1997) show that the fractional content 
of ellipticals is similar to nearby clusters but the SO/spiral ratio is very much lower. This can be naturally 
explained if SO galaxies in rich clusters were spirals in their previous lives before interacting with other 
cluster galaxies and the ICM. When stripped of their interstellar gas by the ICM, SO galaxies could have 
enriched the ICM, perhaps with a strong SNIa abundance signature. If the ICM enrichment were due 
solely to this process, however, the iron abundance in the ICM should decrease systematically with redshift. 
Mushotzky & Loewenstein (1997) find no significant change in the global iron abundance in clusters out to 
redshifts z <~ 0.3, although there is considerable real scatter in the data that may have masked such a trend. 
While the transition of spirals into SO galaxies is appealing as a possible explanation for the large specific 
supernova rates in clusters, morphological classification at large redshift is a difficult art and subject to 
error. In similar studies of distant clusters Stanford et al. (1997) and Andreon (1998) find no depletion of 
SO galaxies at high redshift. 

Another less conventional cluster enrichment possibility is that some of the cluster gas metals has 
come from outflows from dwarf galaxies. At least one author (Trcntham 1994) has claimed that dwarf 
galaxies can supply all of the gas and metal enrichment observed in rich clusters. While few would agree 
with this extreme hypothesis, some fractional contribution from dwarf galaxies to the ICM enrichment 
cannot be entirely dismissed. ^From studies of intergalactic absorption lines in quasar spectra, Cowie & 
Songaila (1998) and Lu et al. (1998) find evidence of an average carbon enrichment of ~ 0.003 solar in 
the intergalactic medium at high redshifts z ~ 3, even in low density voids. Although this value is 100 
times lower than abundances in rich clusters, it does suggest the presence of an alternative source of metal 
enrichment. According to Cowie & Songaila, "early generations of small galaxies might be much more 
efficient at ejecting heavy elements ... than has previously been thought." Lu et al. maintain that the 
variation of the carbon abundance with column (and physical) density in La clouds rules out the possibility 
of (uniform) metal contamination of the intergalactic medium by population III objects. 

In another intriguing recent observation with ASCA, Hattori, et al. (1997) have discovered an X-ray 
cluster at a very high red shift z = 0.92 having luminosity L x — 8 x 10 44 erg s _1 . It is most remarkable that 
this cluster is optically dark. After several sensitive optical searches only the central CD galaxy has been 
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found at this rcdshift. Other cluster member galaxies of average luminosity, if present, should have been 
observed and were not. This cluster is similar in all respects to those at lower redshifts (z 0.4) except 
that its mass to light ratio, <~ 3000, is about ten times that of less distant clusters. Given the apparent 
absence of prominent cluster galaxies, it is most astonishing that the iron abundance in the hot cluster gas 
(kT = 8.6 keV) is zp e — 1.7^ ' 7 4 in solar units. One possible interpretation would be that the cluster gas is 
enriched by a multitude of low mass galaxies which has escaped optical detection. If dwarf galaxies created 
some of the metals in rich clusters, it is likely that this happened at high redshifts to account for the relative 
constancy of the iron abundance found by Mushotzky & Loewenstein (1997) at red shifts z < 0.3. 

Some of these alternative cluster enrichment hypotheses are subject to observational test. We have 
already mentioned that there is a small, but real, scatter in the global iron abundance in cluster gas 
(Mushotzky 1998). Since the SO/spiral ratio also varies among clusters, a weak positive correlation of this 
ratio with gas phase abundance excess would be expected if SO stripping is an important contributor to 
cluster gas enrichment. Abundance inhomogeneities in the hot cluster gas may provide another clue to 
the origin of cluster metals. If only a few bright E galaxies are responsible for most of the intracluster 
gas enrichment, significant abundance inhomogeneities and radial gradients would be expected due to the 
limited previous orbital experience of these galaxies. If the cluster gas was enriched by a larger number 
of dwarf galaxies (most of which may no longer be visible), then abundance inhomogeneities will be much 
reduced. Such studies may become possible as the spatial resolution of X-ray detectors improves. 

The possibility of galactic sources of cluster enrichment other than E and SO galaxies has been 
considerably strengthened by recent ASCA observations of cluster iron and silicon abundances. In a study 
of 40 clusters Fukazawa et al. (1998) find that the ICM silicon abundance increases with cluster richness as 
kT increases from ~ 2 keV to ~ 9 keV, but no appreciable change in the iron abundance is indicated over 
this range. The explanation for this behavior is unclear at present, but it very definitely requires at least 
two sources of ICM enrichment, otherwise the iron to silicon ratio would be the same in all clusters. For 
this range of cluster gas temperatures, kT > 2 keV, the constancy of global gas and iron masses relative to 
cluster Cb suggests closed box environments (Rcnzini 1997). However, the sense of the silicon enrichment is 
opposite to that expected if more spiral galaxies in richer clusters (high kT) convert to SOs. Nevertheless, 
additional sources of cluster gas metals may provide the most satisfactory resolution of the inconsistent 
enrichment histories of massive ellipticals and rich clusters that we have demonstrated here. 

Studies of the evolution of hot gas in elliptical galaxies at UC Santa Cruz is supported by NASA grant 
NAG 5-3060 and an NSF grant for which we are very grateful. In addition FB is supported in part by 
Grant ARS-96-70 from the Agenzia Spaziale Italiana. 
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APPENDIX 

In this Appendix we briefly discuss the evolution of hot gaseous halos in NGC 4472 in the context of 
two different cosmologies: an open universe with f2 = 0.3 and a flat, low density universe with = 0.3 and 
Cl\ = 0.7. When 0,^1 the evolution of the dark matter can no longer be described by the self-similar 
solutions of Bcrtschingcr (1985). Instead, we adopt a simple top hat overdensity perturbation at some early 
time designed to concentrate the dark mass of NGC 4472 by the current time. The dark matter is followed 
as a separate zero pressure fluid which builds a stationary NFW core that grows with time, conserving total 
dark mass. For simplicity we assume the NFW profile for NGC 4472 is identical to that in the flat universe 
described in the main text. 



A. Halo and ISM Evolution in an Q, = 0.3 Universe 

Since all matter in an open universe is unbound, the initial perturbation must be large and non-linear 
in our simple model of galaxy formation. For the Q, = 0.3 universe we introduced a top hat perturbation 
of overdensity (1 + 8)p(ti) with S — 5 extending to Ri = 43.5 kpc at time ti = 10 8 yrs. Baryons with 
mass density SI/, = 0.05 also participate in this top hat perturbation and by time = 1.5 Gyrs enough 
gas is concentrated to form the stars in NGC 4472 in the usual way described in the main text. We 
chose r\n — 0.3r/ st( i for the SNII number (taking e sn — 1) and the Type la supernova rate is specified 
with parameters SNu(i„) — 0.03 and p = 1. In Figure 8 we show with a solid line the gas density and 
temperature distributions after time t n — 13 Gyrs. Although the agreement is not perfect, it is adequate to 
illustrate that the global evolution of extended hot halos in ellipticals in an open universe is quite similar 
to the solution for the flat ft = 1 universe discussed in the main text. 

At 13 Gyrs the transition from the stationary NFW dark halo core to the secondary infall occurs at 
~ 200 kpc, which is just beyond the outermost observations of the X-ray image of NGC 4472. Since the 
potential has a strong slope change at this radius, this feature at the outer dark core might be visible in the 
outer X-ray images of some massive ellipticals if the universe were open. For open universe parameters, i* 
cannot be much larger than ~ 1.5 Gyr otherwise too much cold gas gathers in the perturbation potential. 
In the open cosmology the total baryon fraction is higher at large radii ~ 0.16, similar to rich clusters; 
although the final gas density is necessarily similar to that in the flat universe solution, the dark matter 
fraction is lower. 

We also computed the radial iron and silicon abundances in the gas at time t n = 13 Gyrs and find 
a general agreement with observed abundances, although the fit could be improved with a slightly lower 
rjii. Successful solutions for NGC 4472 in the context of an open universe are incompatible with large 
r]ii rts 2.&rj st d required for global cluster enrichment. 

B. Halo and ISM Evolution in an = 0.3, Oa = 0.7 Universe 

Galaxy evolution in this cosmology was initiated with a top hat perturbation of radius Ri = 3 kpc 
and relative overdensity 5 = 0.59 at time ti — 10 6 yrs. For the supernova rates we take r\n = f] st d and 
SNu(t„) = 0.03 with p = 1. The stars in NGC 4472 are assumed to form at t* = 3 Gyrs. After 13 Gyrs 
the gas density and temperature profiles are shown with dashed lines in Figure 8. The gas density profile 
is slightly too peaked but is otherwise similar to the observed variation. The gas temperature exhibits a 
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curious undulation near r = 90 kpc that is not apparent in the observations. Nevertheless, the computed 
temperature profile for this model agrees with observed temperatures far better than profiles computed 
without secondary infall (see Brighenti & Mathews 1998b). 

Our purpose in illustrating these results with alternate cosmologies is not to achieve the best possible 
fits to the current density and temperature distributions by careful adjustment of the parameters. We only 
wish to demonstrate that such models are feasible and that the principal conclusions reached in this paper 
do not depend on the simple flat cosmology that we adopt. 
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Fig. 1 . — Solutions of Equations (4) - (6) for cluster enrichment parameters in terms of the number of Type 
la supernovae per unit Lb, nj a (bottom x-axis) and the current Type la rate SNu(t„) in SNu units (top 
x-axis). Upper left: Number of Type II supernova per stellar mass rju; Upper right: IMF-averaged iron yield 
for Type II supernovae (yFe.ii)', Lower left: Fraction of metals from Type II supernovae in stars /*. Solid 
lines: correspond to no stellar enrichment by SNIa (g* = 0); Dashed lines correspond to enrichment by half 
of the SNIa ejecta (g* = 0.5). Lower right: Total fraction of iron produced by SNIa F] a (dot-dashed lines) 
and the total fraction of iron in the gas phase produced by SNIa i*/ a ,g {dotted line) when g* = 0.5. Fj a and 
Fi a .g for the 5* = case (solid line). 

Fig. 2. — Upper panel: Variation of the umber of Type II supernova per stellar mass r]u(x) with IMF slope 
x for three pairs of upper and lower mass cutoffs: solid line: (m u , mi) = (100, 0.08); long dashed line: (m u , 
me) = (40, 0.08); short dashed line: (m u , mi) = (40, 0.3). Lower panel: Variation of 0, the fraction of initial 
stellar mass expelled 11 Gyrs after a single burst of star formation. The IMF notation is identical to that in 
the upper panel. 

Fig. 3. — Upper panel: Current gas density profile of the standard model (solid line) and the cluster model 
(dashed line) compared with observations of NGC 4472. Filled circles are NGC 4472 gas densities observed 
with Einstein HRI (Trinchicri, Fabbiano & Canizares 1986) and open circles are densities determined with 
ROSAT HRI and PSPC (Irwin & Sarazin 1996). Lower panel: Current gas temperature profile of the 
standard model (solid line) and the cluster model (dashed line) compared with ROSAT gas temperature 
observations of NGC 4472 (Irwin & Sarazin 1996) shown with errorbars. 

Fig. 4. — Radial variation of iron (upper panel) and silicon (lower panel) abundances for the standard model. 
The ASCA iron observations (from Matsushita 1997) are with the SIS detector (filled circles) and the CIS 
detector (open circles). Solid lines: total abundances; Long dashed lines: abundance contribution from SNII; 
Short dashed lines: abundance contribution from SNIa; Dotted lines: abundance contribution from stellar 
outflow. 

Fig. 5. — Cumulative mass of iron in the hot interstellar gas at time t n = 13 Gyrs as a function of galactic 
radius for the standard model (solid line), the cluster model (dashed line), and the la model (dotted line). 

Fig. 6. Radial variation of iron (upper panel) and silicon (lower panel) abundances for the "cluster" model. 
Notation is identical to that in Figure 4. 

Fig. 7. — Radial variation of iron (upper panel) and silicon (lower panel) abundances for the "la" model. 
Notation is identical to that in Figure 4. 

Fig. 8. — Current gas density (upper panel) and temperature (lower panel) profiles for the open universe 
(il = 0.3) (solid lines) and the low density flat universe (O = 0.3, J7a = 0.7) (dashed lines), plotted against 
observations of NGC 4472 described in Figure 3. 
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